General Disclaimer 


One or more of the Following Statements may affect this Document 


• This document has been reproduced from the best copy furnished by the 
organizational source. It is being released in the interest of making available as 
much information as possible. 


• This document may contain data, which exceeds the sheet parameters. It was 
furnished in this condition by the organizational source and is the best copy 
available. 


• This document may contain tone-on-tone or color graphs, charts and/or pictures, 
which have been reproduced in black and white. 


• This document is paginated as submitted by the original source. 


• Portions of this document are not fully legible due to the historical nature of some 
of the material. However, it is the best reproduction available from the original 
submission. 


Produced by the NASA Center for Aerospace Information (CASI) 





X-732- 73-233 

PRBWNT 

life 70 S??- 

OYNAMICS 

OF FLEXIBLE SPINNING SATELLITES 
WITH RADIAL WIRE ANTENNAS 


R. W. LONGMAN 
J. V. FEDOR 


(NASA-TK-X-70584) DYNAMICS OF FLEXIELE 
SPINNING SATELLITES WITH RADIAL WIRE 
ANTENNAS (NASA) 33 P HC $4.75 CSCL 22B 


SEPTEMBER 1973 


— GODDARD SPACE FLIGHT CENTER 

GREENBELT, MARYLAND 

”' r i: V 



N 7 4- 17572 
Unclas 

G3/31 30799 


X-732-73-233 

Preprint 


DYNAMICS OF FLEXIBLE SPINNING SATELLITES 
WITH RADIAL WIRE ANTENNAS 


R. W. Longman 
J. V. Fedor 


September 1973 


Goddard Space Flight Center 
Greenbelt, Maryland 


CONTENTS 


Pages 

Abstract v 

Introduction 1 

Quadratic Approximation to the Lagrangian 2 

The Linearized Differential Equations 9 

Mode Shapes and Natural Frequencies 11 

Application of Results to IMP-J Spacecraft 16 

Fundamental Spin-Up Equation: Analytic Solution 18 

WKB Approximation 21 

Asymptotic Expansion Solution 23 

Planar Mode Response to Spin-Up 25 

References 28 


ILLUSTRATIONS 

Figure Pages 


1 The IMP-J Spacecraft vi 

2 Inertial and Hub Fixed Coordinate Systems 3 

3 Variables Specifying Position of Wire Element 6 

4 Vibrational Mode Shapes 11 

5 Root Locus Plot 17 


PRECEDING PAGE BLANK NOT FILMED 
Ui 



IV 



DYNAMICS OF FLEXIBLE SPINNING SATELLITES 
WITH RADIAL WIRE ANTENNAS 


Abstract 

A dynamic analysis is presented for a spin stabilized spacecraft employing 
four radial wire antennas with tip masses, a configuration first employed in the 
IMP-J spacecraft. The use of wires in place of the usual booms represents the 
ultimate in weight reduction at the expense of flexibility. The satellite is modelled 
as a 14 degree of freedom system, and the linearized equations of motion are 
found. The lowest order vibrational modes and natural frequencies of the gyro- 
scopically coupled system are then determined. Because the satellite spin rate 
is decreased by antenna deployment, a spin-up maneuver is needed. The response 
of the time varying mode equations during spin-up is found, for the planar modes, 
in terms of Bessel functions and a Struve function of order -1/4. Because tables 
of the latter are not readily available, the particular solution is expressed in 
various forms including an infinite series of Bessel functions and a particularly 
useful asymptotic expansion. An error formula for the latter is derived showing 
that it gives good accuracy. Also, a simple approximation to the complementary 
function is obtained using the WKB method, and the phase error in the approxi- 
mation is shown to be small. 
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DYNAMICS OF FLEXIBLE SPINNING SATELLITES 
WITH RADIAL WIRE ANTENNAS 


Introduction 


A dynamic analysis is presented for a spin stabilized spacecraft employing 
four radial wire antennas with tip masses. Most satellite designs use booms 
which are made as rigid as possible within weight and storage limitations. The 
design analyzed here represents a logical extreme of weight reduction, in which 
the booms are replaced by wires; and correspondingly it represents an extreme 
in flexibility since the wires have essentially no bending stiffness. An effective 
bending stiffness is supplied by the spin of the spacecraft. 

Among the advantages of the wires over the boom configuration are that the 
wires are lighter, and the deployment devices simpler. In a spin stabilized 
spacecraft the deployment decreases the rotation rate. Hence, a spin-up 
maneuver is required. Here again the wires have advantages. Less fuel is 
needed in the spin-up, creating an additional savings in weight, and there are 
no problems of buckling at the root during the maneuver. Prior to planning of 
the IMP-J spacecraft (see Figure 1) for which this study was made, probably the 
major deterrent to the use of the radial wire configuration was the apparent 
complication in predicting dynamic behavior, together with a fear of the effects 
of extreme flexibility on spacecraft behavior and stability. These difficulties 
are more apparent than real for a space vehicle with high spin rate. As is 
evident from the present analysis, the radial wire problem is an order of mag- 
nitude more tractable than the problem with nonzero bending stiffness. 

For purposes of the analysis, the satellite is modelled as a symmetric rigid 
body with four radial spherical physical pendulums. Thus, the model predicts 
the behavior of only the lowest order vibrational modes in which the wires be- 
have as rigid bodies hinged at their attachment points. These modes definitely 
predominate at high spin rates, as in the case of the IMP-J. This vehicle is 
designed to spin at least at 23 rpm at all times. The wires when fully extended 
arc 200 ft. long, and weigh 1/2 gram per foot with small 3 gram tip masses at 
their ends. Note that the distributed mass of the wire is much more important 
than the trip masses, so that a simple pendulum model cannot be substituted for 
the present physical pendulum model with distributed centrifugal loading. 

The quadratic approximation to the Langrangian is obtained, md from it the 
linearized equations of motion for the 14 degree of freedom systen . A 
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trims formation of variables gives these equations in a form where the vibrational 
mode shai>es and their natural frequencies can be recognized. The presence of 
gyroscopic coupling complicates this process, and produces an eighth order 
characteristic equation which must be factored to get some of the system natural 
frequencies. This problem is then reduced to one of finding roots of a cubic 
equation, and an approximate method based on root locus techniques from control 
theory is used to obtain the roots in the case of the 1MP-J. Implicit in this solu- 
tion for the mode shapes and mode frequencies, is the general solution of the 
differential equations of motion oi the satellite for any given initial conditions. 

The spin-up maneuver needed in antenna deployment is a planar operation. 

In order to analyze the behavior of the wires during spin-up, the differential 
equations for the planar mode variables are obtained, linearized about the rigid 
body spin-up solution. Although these equations are linear, they have time 
varying coefficients, and thus present some difficulties in obtaining their general 
solutions. By a transformation into the complex plane, the gyroscopically coupled 
modes, as well as the other planar modes, are all shown to be governed by one 
fundamental equation. By proper transformations of both the dependent and 
independent variables the analytic solution of this equation is obtained in terms 
of Bessel functions of order ±1/4 and a Struve function of order -1/4. 

Because tables for the needed Struve function are not readily available, ap- 
proximate solutions are obtained which give additional insight into the satellite 
behavior. The complementary function can be approximated using the WKB 
method, and the phase error is shown to remain small. The Struve function can 
be approximated using a Taylor series, but convergence is very slow. Some- 
what faster convergence is obtained by deriving an expansion in terms of an in- 
finite series of Bessel functions. However, it is shown that the most useful 
expression for the particular solution is in the form of an infinite asymptotic 
expansion which diverges for all finite time, but nevertheless gives very good 
accuracy when a small number of terms are used. An error formula is derived 
which expresses the error present in this approximation, and also gives a lower 
bound on the number of terms which can be included before the error begins to 
increase. Combining these approximate solutions gives the general solution for 
the planar response of the system under a spin-up torque. 


Quadratic Approximation to the Lagrangian 


We v i h to obtain the equations of motion linearized about a steady rigid 
body spiL for the modal analysis, and linearized about a rigid body spin-up 
oor e?po g to a constant applied torque for the spin up analysis. The formu 
lation usvd will obtain both sets of linearised equations simultaneously. It 
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Figure 2. Inertial and Hub Fixed 
Coordinate Systems 


suffices to use the quadratic approximation to the Lagrangian in order to obtain 
linearized equations of motion using Lagrange’s equations. Since there is no 
potential energy we need only the quadratic approximation to the kinetic 
energy, T. 

Figure 2 shows an inertially fixed coordinate system XYZ and coordinates 
xyz fixed in the symmetric hub of the spacecraft and centered at its center of 
mass. Let R ov be the vector from the center of the inertially fixed coordinates 
to an arbitrary volume element dV of the spacecraft, and let R 0H and R HV be as 
shown. Then the velocity of the element dV relative to inertial space is 
dRov/dt ,, the time derivative in inertial (I) axes. The spacecraft kinetic 
energy can be written 



where p is the density and m t is the mass of a tip mass. 
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The angular velocity <*> of the hub axes relative to inertial space expressed 
in hub coordinates will be needed. It can be written as 


W 1 


• • 

9 cos 0 ^ cos 0 3 + 0 ^ s i n 6 3 

^2 

- 

- 9 cos 0 2 s i n 9 3 + 9 ^ cos 0 3 



V in9 j*^ 


where , , 9 , and ; are 1, 2, 3 type Euler angles, i.e. a rotation through 
angle ' ' in a right handed sense about X, followed by a rotation d 2 about the 
resulting Y axis, and then a rotation 6 about the resulting Z axis gives a set 
of axes parallel to the hub xyz coordinate system. 

To perform the linearization we can write 


; 3 = V*) + 

& 3 - \p( t) + be 3 


where and 8^ are assumed small quantities. For the modal analysis we 
linearize about the steady rotation of a rigid body, so that w 0 (t) = a constant 
angular velocity about the spin axis, and i/»(t) = ^ t. For the spin-up analysis we 
need to linearize about the values of 9 and 9 3 corresponding to a constant 
torque Q applied about the spin axis, assuming the satellite to be rigid and to 
have moment of inertia I about that axis. Thus, we let 


vn = % + (Q/V 1 
Wo = V + i (Q'lj ) {2 

where is the initial spin rate before spin-up. 

The quadratic approximation to the angular velocity can now be written 
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4 • • • 


0 



COS 4 (' 

sin i p 

0 

+ 


s i n 4' 4 IL 
T 2 

cos \p 

A 



A 



W n 4 W, 4 w 2 4 . . . 


- 0 si nf 4 $ i >^ 3 cos 0 
~ ^1 ^^3 cos 0 ** ^2 s ‘ 


Here and elsewhere the three dots indicate terms of higher order than quadratic. 
The other angles and , are expanded about zero nominal values. 

The true generalized coordinates for the center of mass of the hub will be 
the components X, Y, Z of R OH in inertial space. Their derivatives X, Y, Z 
represent the components in coordinate system I of the inertial derivative 
dU OM /dt I . It will be convenient to define x,,, y H , z H as the components of this 
vector in the hub H coordinate system (not to be confused with components of 
the derivative of H ol| as seen in H coordinates, dR OH /dt H ). 

The first term in Kq. (1) represents the kinetic energy of the hub, and can 
be written as the sum of the translational kinetic energy* of the center of mass 
and the rotational kinetic energy about the enter of mass: 




+ + 


1 I UJ 2 

2 HJ 3 


where m H is the mass of the hub, and I H1 , I |U , I H3 are the principle moments 
of inertia about x, y, and z axes respectively. T v e z axis is the axis of sym- 
metry. The quadratic approximation for this is obtained by substituting the 
above quadratic approximation to w. 

In order to calculate the second and third terms of Eq. (1) we must de- 
termine the ire''tial velocity dH ov /dt , of a volume element of the wire or the 
tip mass. Using Coriolis law we can write 


dR OV 

d ®OH 

^-HV 

. dR o„ 

(, ^HV 

dt 

(It 

I 

dt 

I 

dt 

1 

dt 

I 



The term dR HV /dt | represents the rate of change of R HV as seen in H 
coordinates. 

Figure 3 shows a volume element dV of the i ,h wire. The variables a, and 
b, are the in-plane and out-of-plane angles specifying the deviation of the wire 
from the nominal position. The c ( represent the angle from the x hub axis to 
the i ,h nominal wire position, and c | is taken as zero. From the figure, the 
x,y,z components of the volume element position are easily written. Then the 
quadratic approximations to the matrix [ R HV ] H of components of in H 
coordinates can be written 



cos c 

1 


- n ( sin c ( 


cos c 

1 

I r „ v ) h 4 r ) 

sine. 

i 

+ q 

a. cos c. 

■ i 

- \ + a ?) 

sin c . 

t 


0 


b 

_ i _ 


0 


U 0 ♦ U, ♦ u 2 


+ . . . 


i 



Figure 3. Variables Specifying Position 
of Wire Element 
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Similarly, the matrix [dR HV /dt| H ) H °* derivatives of x, y, z components is 



- n. sine. 


COS C t 

dR HV /d, l H ) H 5 <> 

a. cos c. 

i i 

- q(b.b. 4 a,n. ) 

sin c i 


b. 

i 


0 


U, + U 2 ♦ 


Combining these results, the matrix Idl^/dt , l H of components in H coordinates 
can be written 


ld B OV /dt ll = ^H Z H ]T + + d 2 + fW 0 4 W, 4 W 2 ) (U 0 


♦ u. 


U 2 )4 


where superscript T indicates transpose, and the tilde indicates the matrix 
equivalent of a cross product (if the components of Wj are w ( , then the i, k 
component of the square matrix \fr, is given by 


3 

f ijk w j where jk = y(» - j) ( j - k) (k - i)). 
j' i 


The quadratic approximation to the kinetic energy, per unit mass, of the 
volume element is obtained by calculating the dot product 


5 l dR ov 


di j] H [dR oy dt j H 


and retaining only terms through second order. Multiplication by m t and 
substituting i for q (see Fig. 2) gives the kinetic energy of tue i th tip mass. 
Mu)*'plication by the density per unit length of the wire, p , and integrating q 
from 0 to i gives the kinetic energy of the i th wire. Adding these results for 
each wire and tip mass t^ the kinetic energy of the hub gives the 
quadratic approximation of the total kinetic energy 
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T = 3 M (*« * *1 * *3) ♦ W*% 4 


4 3 I i<‘ : i 4 ' ! j> 4 \ V s ®*' 2 4 M ' <w « 4 s ®*> 2Z ** 

i * 1 


t l M 3 W »£] a ? 4M 2 i H T}i 
• * 1 1*1 

4 4 

- ^ M i w o ^ (a ? + b i> + 2 “ 3 jZZ“‘ + ^ ) 
i * : i*i 


♦ M 2 l- x H ( a. sinr + w Q a. cos c. ) + y H ( ir cos c. - 


> r l 


4 

4 M 1 Z^ ^1 sin ('p + c i) - w o b i cos(0 + c. )] 

i » 1 


- " 2 [b i e os(</> + c^ 4 w 0 b. sin(</> + c. )] } ♦ . 


where 


M = «^ + 4m + 4m t 



w Q a. sin c.)] 
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M 


2 



M, = M 3 + rM 2 


*1 “ 1 + 
*3 = *H3 + 


4 ^m t + 1 mj -t 2 + 8 ^m t + 1 m j r-t + 4(m t + m) r 2 
2 ^m t + ^ mj 't 2 + 4 ^m t + ^ mj r'i + 2 ('m t + m) r 2 


Ip Ij , and 1 3 are the three principle moments of inertia of the entire spacecraft 
about the center of mass assuming it is rigid; I H i , Ih 2 » I H 3 are the correspond- 
ing inertias for the hub; r and i are defined in Figure 2; m is the mass of a 
wire, m t that of a tip mass, and that of the hub. 


The Linearized Differential Equations 


The linearized differential equations are obtained by applying Lagranges 
equations to the above kinetic energy for each of the 14 generalized coordinates 
X, Y, Z, 0j, 0 , jf? 3 , a, , bi . Nqte that X, Y, Z are cyclic so that the equations 
reduce to dT/BX = Cj, 3T/3 Y = C 2 , BT/3 Z = C 3 , C k all constant. Further- 
more, we can chose the I axes so that there is no momentum of the spacecraft 
relative to I at time t = 0. Assuming this has been done, the constant generalized 
momenta are zero, Cj = C 2 = C 3 = 0. 

Direct application of Lagranges equations yields linear differential equations 
with variable coefficients. Some of this difficulty can be circumvented by keep- 
ing the Xn, y H , and z H coordinates. We recognize that to a linear approximation 


BT _ BT BX + BT BY + BT BZ 

Bx rf 3X Bi^ BY BXjj BZ B^ 

BT , BT . . 

= — cos^ r s i n i p + . . . 

dX BY 

or BT/Bx h = 0, and similarly for y H and z H . 


9 


• * 7 ~ 


The linearized equations of motion are then 


4 

- M 2 (a. sin c. + w Q (t) a. cos c £ ) = 0 

i = l 


(2a) 


r H T ‘"2 C ° S C * " 3 i Sin C i> = 0 

i=l 


My u + M, 


(2b) 


Mi « + M 2 ^ b. = 0 


(2c) 


4 

! 1^1 + l 3 + V 1 )^ +M i^Ubi + w2b £ ) sin(>//+ -w 0 b. cos^+cj = Q(2d) 

i=l 


1 1^2 ~ J 3 W 0 ( * ) " M 1 / . + w J b i) cos(0 + c.) + w 0 b. s in(0 + c. )] =0 


(2e) 


i=i 


i 

V*, + M i ) “i = o 

i=l 


( 2 £) 


M 3 a i * M 1^3 + rM 2 w o a i *V(-VVh> sine, *(y s t» l) i,)cosc 1 l =-M,w 0 (2g) 


M 3 ii, +Mj'i tM^bj 3M 1 [(e i + 2w o ^)sin(^tc.)t(-9 2 + 2w o 0 l )cos('/- + c i )] =0 (2h) 


where i = 1, 2, 3, 4. 
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Mode Shapes and Natural Frequencies 


We wish to find a change of variables which will uncouple the equations of 
motion, and thus determine the normal modes of vibration and their natural 
frequencies. Because of the spin of the spacecraft, we cannot actually expect 
complete decoupling. Certain modes will be gyros copically coupled. It might be 
possible to determine the needed mode variables mathematically, but it is much 
easier to appeal to physical intuition to guide us. Figure 4 presents eight ways 
in which one would expect the system to oscillate. Associated with each mode 
shape is a mode variable which should behave sinusoidally 



Figure 4. Vibrational Mode Shapes 
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“l = \ <*1 * a 2 + 8 3 + V 

°2 = * S 2 + °3 - V 

“3 = \ < a i + a 2 - a 3 - 3 4> 

“« = \ <- a i + a 2 * a 3 - a «> 


^1 = ^ b l + b 2 * b 3 + b 4> 

V, = \ < b l - b 2 + b 3 - b 4> 

^3 = | < b l + b 2 ' b 3 - b 4> 

^4 = 5(- b l ,b 2* b i- b «> 


Note that a 2 and f3 should be completely Independent of other mode variables. 
Also, cij and /3 equations will couple with the rotation and translation equations, 
but will be independent of other mode variables. However, we expect a 3 and 
a 4 to be gyros copically coupled, and fi 3 , to be coupled with the precession 
and nutation of the hub. 

Setting w 0 = co 0 , 4> = o. 0 t, c k = (i - l)rr/2, and using mode variables, Eqs. 
(2a-c) and (2f) become 


x,, = (2Mj/M) [(a 3 + i 4 ) + o> 0 (a 3 - a 4 )l 
y H = - (2Mj/M) [(i 3 - a 4 ) - ^ 0 (a 3 t a 4 )] 


i„ = - (4Mj/M) 


8^ = - (4M, /Ij) a, 


Adding Eqs. (2g) for [ = 1, 2, 3, 4 and using the above S3 , and similarly adding 
Eqs. (2h) and using z H gives 

[M 3 - (4Mj/I 3 )] Oj + rMjwJaj = 0 

[M 3 - (4Mj/M)] \ = 0 
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Therefore, the natural frequencies associated with these modes, normalized by 
the rotation rate , are 


Q al = {rM 2 /[M 3 - (4Mj/I 3 )]} 1/2 


fig, = {Mj / [M 3 - (4M 2 /M)] } 1/2 


Similarly, taking the linear combinations of Eqs. (2g) and (2h) indicated by the 
definitions of a 2 and (3 give 


n a2 = {rM 2 /M 3 } 1/2 


and 


n 


02 


{M/M 3 } 1/2 


n 2 

a 2 


cv 2 cl 

0 2 


= 0, 


6 + n 2 co 2 3 = o 

2 02 0 2 


Adding Eqs. (2g) as indicated by a and a , and using x , y , gives two 

3 4 H n 

equations 


[M, - (2M? /M)] a. 


[4M 2/ M] 


a 

0 4 


+ [rM, + (2M 2 /M)] 


oj 2 a 
0 3 


= 0 


[M 3 - (2M 2 /M)] a 4 - [4M 2 /M] co Q a 3 + [rM 2 + (2M 2 /M)] 


U 0 3 4 


= 0 


with skew-symmetric gyroscopic coupling of the first derivatives. The normal 
ized frequencies associated with the roots of the characteristic equation for 
this system are 
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where 


n 


a 3 . 04 



+ 


C ) 


1/2 


B = 4M2/(MM 3 - 2M j ) 


C = (rMM 2 + 2M 2 )/(MM 3 - 2M*) 


To uncouple the equations for (3 and /3 4 , obtained using Eqs. (2h), from 
the remaining equations we must use Eqs. (2d) and (2e). Setting w 0 = u Qt 
0 = w t, Cj = (i - 1 ) 77 / 2 , and using fi 3 and /3 4 to eliminate the b i , the latter 
equations still contain time varying coefficients which complicate the elimination 
of and from Eqs. (2h). However, we note that when w Q = oo Q these two 
angles are cyclic coordinates. We then define a rotation of the derivative 
variables which will eliminate the explicit time dependence. Let 


1 7 1 = ^ cos (%t) + ^2 sin(w Q t) 


V 2 = " $2 + ^2 COS (% t > 


and Eqs. (2d) and (2e) become 


I 1 ^i + ( l 3 ~ V ^0^2 + 2M i [ ^3 + ^4) + "0^3 + = 0 


V2 " ( J 3 " x i> V^i ~ 2M itC/3 3 - \) + wJ(/3 3 - /3 4 )] = 0 


(Note that these same equations might have been obtained more directly by 
using quasicoordinates ojj, a j , cu 3 from the beginning.) The equations for Z 3 
and ? 4 from Eqs. (2h) in terms of rj l and t] 2 are 
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M 3' 3 + + V + 4 V ] = 0 

M 3 4 4 M 1^0^4 + \ M 1 trt, + V ‘ M ^ 0 


The natural frequencies associated with J3 and (3 must be determined from the 
characteristic equation of the above system of four simultaneous equations: 


{(s 2 + Pa> 2 ) (s 2 + Kai 2 ) - 2PN(s 2 + oj 2 ) (s 2 + K^ 2 )} 2 


+ {2PNo; 0 s(s 2 + a.' 2 ) (K - l)} 2 = 0 


where 


P = Mj/M 3 . N=M,/Ij. and K = (I 3 - I l )/I 1 - 


It would appear that we must solve for the roots of an eighth order polynomial. 
The normalized frequencies , 4 are the positive roots Q of this equation 
when we set s = However, we note that the equation is then of the form 

x 2 - y 2 = (x - y) (x + y) so we can reduce the problem to finding the roots of two 
fourth order equations 


(H 2 - P) (P 2 - K 2 ) - 2PN(P 2 - 1) (P 2 - K) 


± 2PNQ(P 2 - 1) (K - 1) = 0 


However, from the form of these equations, if 0* is a root of one, then - 0* is a 
root of the other. Therefore, we need only solve for the four roots of one of the 
quartics and take their absolute values. We can still do better if we observe that 
the precession frequency (or its negative), Q = |K , is a root of the quartic. We 
conclude that the four frequencies fi- 3 4 associated with and /? 4 mode 
variables are 
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V/w = lla- 1 .!". 


and the absolute values of the three roots of 


(Q 2 - 1) [fi - (K - 2PN)/(1 - 2PN)J - k(ft + K) = 0 


where k = (P - 1)/(1 - 2PN). Approximate values for these three roots can be 
obtained easily using root locus techniques as shown below. 


Application of Results to IMP-J Spacecraft 


The parameters of the IMP-J spacecraft are: M t = 3 gms, m = 100 gms, 
i = 200 ft, r = 2 ft, M = 19.37 slugs, I H3 = 117 slug ft 2 (hard booms deployed), 
and the transverse inertias are approximated by 1„, = 76 slug ft 2 . Observe that 

H 1 


P = 


1 + (rAf) 





and is near unity (actually 1.01459) because of the ratio (v/i). The factor k is 
then small (0.0539). The above cubic equation can be written as 


x _ k(fl + K) 

(fi 2 - 1) (fi - G) 

where G = (K - 2PN)/(1 - 2PN). It is thus in the form needed for a root locus 
plot; there is one zero at -K, and three poles at ±1, and G, and we are 
interested in the roots for a small negative gain (-k). For the IMP-J K = 0.875, 
G = 0.540, and the root locus plot is shown in Figure 5. An approximate value 
for each root can be obtained by calculating the root sensitivity at k = 0, i.e. 
finding the rate of change of the root location with k at k = 0, and extrapolating 
along the tangent: 


fi(k) =fl(0) + 



k 
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-1 - 0.87 | + 0.54 
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♦ 1 
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— ► 


Figure 5. Hoot Locus Plot 


where 


dft ft + K 

dk 3 ft 2 - 2Gft - 1 - k 


As a check on the result, the value of k corresponding to the approximate root 
can be calculated by substitution in the root locus equation after solving it for 
k. If the value is not as close as desired, iterating by making a small change 
in ft and recalculating k quickly gives the needed result. 

Using this method to obtain the frequencies for fl /33>/ g 4 and the formulas of 
the previous sect. >n for the other modes gives the following natural frequencies 
normalized by the spin rate : 


ft al = 0.2563 
ft a2 = 0.1208 

ft„, . = 0.1236 and 0.1225 
ft^j = 1.007819 
ft^ 2 = 1.007267 

ft , 3 = 0.8754. 1.002215. 1.09185, and 0.4499. 
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For a spin of 23 rpm the period for the cij mode is roughly 10 sec., while the 
,! 2 , a -j and a 4 modes all give approximately 21 sec. The out of plane modes give 
periods between 2 and 3 sec. with one exception of 5.8 sec. associated with one of 
the / 3 .^4 frequencies. 


Fundamental Spin-Up Equation: Analytic Solution 


We now turn our attention to the dynamic behavior of the satellite during 
the spin-up operation. From physical considerations it is clear that the mode 
« ( is the most directly affected by the spin-up torque. We therefore examine 
this mode first. 

The differential equation for a l is obtained by taking one fourth the sum of 
the four Eqs. (2g) and eliminating 50 using Eq. (2f). The w 0 (t) is taken as 
a o + (Q/I 3 )t. The result is 


a l + ^al 2 '0 + € t) 2 a | - - € ( kjjjMj / r 


where t = Q/(I 3 *- 0 ). The spin-up torque for the IMP-J is 0.96 ft. lbs., which is 
small compared to the total vehicle inertia I 3 , and therefore € is a small 
number. The fi fll is the natural frequency normalized by the spin rate before 
the start of the spin-up, u. 0 . The differential equation can be written in an alter- 
nate form by defining a new time variable r and letting prime indicate differ- 
entiation with respect to r 


r = 1 + e t 


a*; + n „, 2 r’a, = -« [M/HrMj)] (3) 


where fl al - 0 al v« (we will use analogous definitions for all other modes). 
Despite the simple appearance of this equation, it is difficult to obtain the solu 
tion in a form which is easily applied. 

Using the change of variables from a p T to J, £ where 
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£ = - 0 . r : 

^ 2 a 1 


J(f) = (r)- l/J a, (r) 


the differential equation is transformed to 




d£ 


dc 




lM,/(^rM 2 )] C ^) 3 4 


The complementary function for this equation is a linear combination of Bessel 
functions of order ±1/4, and hence the complementary function for is 


a 


lr 


C 


vr 


J. 


4 




( 4 ) 


(c, and c 2 arbitrary constants). Bessel function tables for these fractional 
orders can be found in [1]. 


To obtain the particular solution, note that the Struve function of order v , 
H (£), is a particular solution of 


£ 2 


llJ 

df 2 



(£ 2 -* 2 )J 



(see !2J). By making the appropriate adjustment of constants the particular 
solution for a j can be written as 


“,p = -« 


5 4 




( 5 ) 
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Tables of Struve functions of the needed order are not readily available. One 
can obtain an expression for H. in terms of a product of a fractional power of 
- times an infinite series in <f by using the method of Froebenius on the differen- 
tial equation for H_, 4 ,121. The resulting particular solution for a is 



which represents an entire function. 

Note that 0 al is a large-number, and that the values of t of interest are 
greater than 1. Hence, although the above series coverges for all r, the number 
of terms needed to calculate a good approximation to a is quite large. A 
scries with improved convergence properties can be formed using a Bessel 
function expansion. For any value of n not equal to a negative integer or zero 
12 ] 


(I . y 1 (n + 2k) ry t k) 


Then the driving term to the Bessel equation form of our differential equation 
can be rewritten giving 



Substitute a solution of the form J = 1 a and observe that the 

left hand side of the equation equals [^ 2 - 1/16J J when is substituted for 
J. The coefficients a n are obtained immediately by equating coefficients of 
Bessel functions of like orders 
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This expression represents the exact solution and has better convergence prop- 
erties than the power series solution. However, we will find that for most pur- 
poses an approximate solution in the form of an asymptotic expansion is most 
useful. 


WKB Approximation 

The complementary fu"~iion a found above appears not only in the 
mode, but in other modes as well. For purposes of making rapid calculations, 
as well as to give a somewhat better intuitive feel for the mode response, an 
approximate solution is generated by the WKB method 3i. 

The WKB functions 


W f (T) r [ f (r)J 4 



‘ f (t/) d.*/ 


are easily shown to be exact solutions to 


W; + (f(7) + g(r» W 4 = 0 


1 f" 

g( T ) = - — 

V ; 4 f 



and therefore approximate solutions to aj + f(r) 0 provided is small 
conpared to f( ~ ). For the present problem we can take 


( 1 /r) exp 


ii-fi 

L 2 
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and the complementary function can then be written as 


a , c ( T ) ' A ^ cos (j n *i T2 4 8 ) 


( 8 ) 


wnere A and are arbitrary constants. The function g(r) = -3/(4 r 2 ) is small 
compared to f(r) = fi^r 2 because of the presence of a factor l/e in . 

Applying the well known asymptotic form (for large r) of the Bessel lunc- 
tions would give the same result, and therefore accuracy is good for prcoerly 
chosen & when r is large. The error in this WKB solution when it is carried 
in to small values of r can be investigated by writting 


a i c = C + (r) W + (r) + C.(r) W_(r) 


where the variation of the constants C t = (A/2)e tl accounts for the difference 
between the true solution and the WKB approximation. The rate of change of C t 
with r is given by ( 3j 


uc, ■>: _ 

-j— = 1 — ^ C t 4 C? ex P[+ i ^ al 7 ' 2 ]} 

dT 8fi al r3 


II the C t do not change much, we can calculate the change AC t by integrating 


AC 


t 



±?i * — > 

{1 + exp [7 ifi r 2 7 2 i 5] > dr 

1617 ,r‘ 

1 


treating is a constant. The second term in the brackets is bounded by one 
in magnitude, so that we obtain an approximate expression 


lACyC,; < 3/(8ff„) = 3«/(8^ 0 n al ) 
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. .'f 


showing that the relative error is small even for r as small as T = 1. Thus, no 
significant phase error in & accumulates in our approximate solution even after 
an arbitrarily large number of wavelengths. 


Asymptotic Expansion Solution 


An approximation to the particular solution a lp can be obtained in the form 
of an asymptotic expansion. We formally substitute 


CD 



n - 0 


into the differential equation and solve for the coefficients c n . By the ratio test, 
the resulting series is seen to diverge for all finite r . However, for properly 
chosen N, the series truncated to N terms forms an asymptotic expansion which 
will be shown to give good accuracy 


eM, 


a ~ — 

IP 


% r !V 2 


n- l 


N (- 1 )" 7T (4k + 2 > < 4k + 3 > 

"L- ' 


k = 0 


n = 1 


Q 2n r 4n 

al 


( 9 ) 


Denote the right hand side of this equation by a 1NI and determine the differential 
equation for which it is a particular solution by calculating a" N + Q 2 t 2 a ^ . 

It is seen that the forcing function differs from the desired forcing fimctionby 
the addition of a term depending on r“ 4fN+1 ) . Then the differential equation 
satisfied by the error, Aa = a )P - a JN , is found by subtracting the equation for 
a JN from Eq. (3) 

N 

\ TT (4k + 2)(4k + 3) 

Aa' 4 ? 5 » r»Aa = (-1)"(_L ] 

WV fj2N T 4N44 

\ / al 


The error in the asymptotic expansion is that particular solution Aa satisfying 
Aa = Aa' = 0atT =1. Denote the driving function in the above equation by F(r ). 
Then the desired particular solution is found by variation of parameters to be 
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A a = 



tei(e) <V T ) - ‘V'O <*> 2 (f)J tF(0/W(c)] 


where 0, and <£ 2 are the two linearly independent solutions to the homogeneous 
equation 



and W is the Wronskian. Then^J + x = 0, and <£" + = 0; 

multiplying the first by 0 2 , the second by 0 , and subtracting gives an expres 
sion which when integrated shows that the Wronskian is constant. Let 
1/2 flj r 2 = 1, and use the fact that the Wronskian of J I/ (z), «!_„ (z) is -2 sin 
(un)/(n z ) from [2J, to determine the constant 


*’V[i w (jv) t-i/ 4 (| n a. T2 ) - j;.-. {\ <V J ) .r. I/4 (in al rjj 


2[r, /4 (i) i.' 1/4 (i) - j; /4 (i) .i. 1/4 d)] 


= 2 3 j n 


The Bessel functions are bounded by [2] 


I.W X >I * 1 I.T-,/« (x >l * l/ liH* 


Therefore, the error in the asymptotic expansion solution is bounded by 
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|Aa| <; J [|0 1 (^)| |0 2 (-)| + |0j (r)| k 2 (Oll C|F(c )| /W) d£ 


r 


/ 1 \ 

( SR 

J,/4 

J-. 4 (j n «. T2 ) 



/ 


'.M (j^) 


[|F(^)| /w] d<? 


N 


e 2N+ *• 25 "M, (4k + 2) (4k + 3) 


t | 1 ' ' 

k = 0 

t 

i 

[i 1 

2a.»*>.25 rM2n »I.O 2S p j 


i i 

4N + 2.5 

T 4N ♦ 2.5 


/T 

+ 

i . _i 

] 

4N + 3 

r 4N ♦ 3 
L ' J 

J 


Because e is quite small in the present application it is found that using 
three terms in the asymptotic expansion gives very good accuracy. 


Planar Mode Response to Spin-Up 


Because the spin-up maneuver is a planar one, we now obtain the dynamic 
response of the remaining planar modes under the spin-up torque. Taking the 
appropriate linear combination of Eqs. (2g), letting w 0 (t) = % + (Q/I 3 )t, and 
changing from t to r as the independent variable gives 


n 2 r 3 

a 2 


a 2 = ° 


and the corresponding solution is 


a 


2 
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where fi a2 = Cl a2 co Q /e . The differential equations for the coupled modes a 3 and 
a 4 are found by taking the appropriate linear combinations of Eqs. (2g), using 
Eqs. (2a) and (2b) with a = (i - 1) */ 2 and the above w Q (t), and changing from 

ttO r 


a" + Bra' + Cr 2 a + -8a = 0 
j ^ j 2 ^ 

a" - Bra' + Cr 2 a - —8a = 0 

4 3 4 2 3 

where 


8 = 


4M 2 a> 


o 


e(MM 3 - 2M 2 ) 


p _ %C r MM 2 + 2M 2 2 ) 
e 2 (MM 3 - 2M 2 ) 


Note that this time there is not only gyroscopic coupling of the first derivatives, 
but coupling in the undifferentiated terms as well. In order to obtain an analytic 
solution to this system of equations, use a complex valued change of variables 
to write the system as a single equation 


a 


n 


i8ra' + 



j iSj- a = 0 


where a = a 3 + ia 4> A second change of variables a(r) = v(r)p(r), where 


p(r) = exp 



is chosen to eliminate the first derivative term in the differential equation, 
gives 


which is the fundamental equation which has already been solved. Then the solu- 
tions for a 3 and a 4 , in terms of real valued arbitrary constants, are both of the 
form 



A simple calculation shows that the WKB approximation is 


a 3- a 4 = jz cos ^ O a3 r2 + tjj + -1 COS^I Q m t 2 t S 2 j 


where A Jf A 2 , 8, , and & 2 are arbitrary constants determined by initial 
conditions. 

Although this section is concerned with the planar modes, it is interesting 
to note that two of the out of plane modes are governed by the same fundamental 
spin-up equation. By a derivation closely paralleling the deviation of the mode 
frequencies for / 3 ^ and /? , the spin-up equations are found to be 
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h = c ,,/Tj 


1/4 
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